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ABSTRACT 

The NASA Kepler mission has revolutionised time-domain astronomy and 
has massively expanded the number of known extrasolar planets. However, 
the effect of wide multiplicity on exoplanet occurrence has not been tested 
with this dataset. We present a sample of 401 wide multiple systems con¬ 
taining at least one Kepler target star. Our method uses Pan-STARRS 1 and 
archival data to produce an accurate proper motion catalogue of the Kepler 
field. Combined with Pan-STARRS 1 SED fits and archival proper motions 
for bright stars, we use a newly developed probabilistic algorithm to identify 
likely wide binary pairs which are not chance associations. As by-products 
of this we present stellar SED templates in the Pan-STARRS 1 photometric 
system and conversions from this system to Kepler magnitudes. We find that 
Kepler target stars in our binary sample with separations above 6" are no 
more or less likely to be identified as confirmed or candidate planet hosts 
than a weighted comparison sample of Kepler targets of similar brightness 
and spectral type. Therefore we find no evidence that binaries with projected 
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separations greater than 3,000 AU affect the occurrence rate of planets with 
P<300days around FGK stars. 

Key words: binaries: visual, astrometry: proper motions, stars: planetary 
systems 


1 INTRODUCTION 


'he Keyler 


flRowe et al. 


2015 


Uissio n fiBorucki et al.l 12010}) has to date identified 3697 exoplanet candidates 


and references therein). The confirmed candidates and inferred false pos¬ 


itive rates have been used to determine the faction of stars with shorter period planets 

-1 

3|) with particular a 


(iHoward et ah 


2011 


Fressin et ah 


2011 


Earth-like planets around M dwarfs flMorton &: Swift 


2014 


dention paid to the number o f 


Dressing &: Charbonneaull2015[) . 


The determination of the parameters of transits exoplanets dependents heavily on the ob¬ 


served properties of the host star. Hence the study of exoplane 


and temperature has been a fruitful endea.v our for several groups (iMann et ahl 


hosts’ metallicity, radi us 


2012i 


Muirhead et al. 

2014; 

Newton et al. 

2015) 


2013 


a 


20151) . One key problem with late-type stars is the sig¬ 


nificant pollution of the dwarf s equence by backgr ound giants. While these can be readily 
identified through spectroscopy (iMann et ahl 120121) this is hard to do for large numbers of 
stars. Because planets are difficult to impossible to find around giant stars using transit 


searches fane 


Bastien et ah 


because evolved stars can cause errors in planet property determinations, 


20141 IWang et ahl 120151) , plane t occurrence calcu lations require an estimate 


of the number of evolved stars in the survey fiMann et al 


20121) ■ Therefore, calculating the 


planet fraction requires an accurat e characterisation o f the tens of thousands of Kepler tar¬ 


get stars with no detected planet. 


Mann et al 


fj2013al) showed that reduced proper motion 


is an efficient dwarf/giant discriminator for late K and M dwarfs which can be used on much 
larger samples of stars. 

The quality of proper motions in the Kepler field varies with the magnitude of the target 
star. For the brightest objects either Hipparcos (to U 9) and Tycho (to V ~ 11.5) provide 
proper motion measurements with uncertainties of the order of a few millarcseconds per 
year. Astrographic observations in the UCAC 4 catalogue provide proper motion errors 
below 10 milliarcseconds per year for stars down to i? ~ 16. This covers the majority of the 
Kepler targets but will miss any faint common proper motion companions to these Kepler 
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Wide binarity and planet occurrence in the Kepler field 


targets. Digiti sed archival p 


and USNO-B ( Monet et al. 


lotom etric surveys such as SuperCOSMOS flHamblv et al.ll200ll) 


20031) provide proper motions for these fainter stars. While these 


surveys have the advantage of long time baselines, their photographic nature means that each 
individual epoch is of poorer astrometric quality than modern CCD surveys. Additionally 
the most extensive proper motion survey based on photographic plates, the Superblink 
survey (with errors of approximately 8 milliarcseconds per year) only has a few subsets of 
ata which are public ly available: high pro per motion stars in the Northern hemisphere 


flLepine fc Shara 


2005h and bright M dwarfs ( Lepine &: GaidosIPOllI) . 


Wide stellar companions are common in the field, wit 


having co mpanions wider th an 100 AU flRaghavan et al 


1 >25 % of nearby solar type stars 
20101). Among young stars, disk 


frequency (iKraus et ahl 120121) and disk mass (jPascucci et al 


20071: 


Harris et al. 


2012) are 


suppressed for close binaries but the components of systems wider than 100 AU have similar 
disk properties to single stars. This s uggests that wid e bina ries play a minor role in the 
initial sculpting of planetary systems. iRaghavan et al.l (1200611 found 24 companions wider 


than 100 AU amongst a sample of 131 radial velocity exoplanet hosts. This equates to a 


m 


Raghavan et al. 

(2010 

) which have companions wider than 100 AU. Finally, 

Wang et al. 


(120141 ) studied binary systems closer than 1500 AU hnding a significant deficit of planets 


around the closest binaries with the effect dropping significantly becoming very small for 
systems wider tha n 100 AU. 

Recent work by 


Kaibetah 


(120131 ) suggests a different mechanism by which binaries with 


semi major axes of 1,000 AU or so can affect planetary systems. While orbiting through the 
Milky Way’s gravitational potential, the orbital parameters of wide binary companions can 
be modified by gravitational interactions with passing stars driving them to more elliptical 
orbits and sending the companions close to their primary at periastron. Such close passage 
can disrupt the outer planets in such systems and induce planetary migration and/or ejec¬ 
tion. The magnitude of this effect can be tested by comparing the planet occurrence rate 
around wide binary components with that around field stars. However this is complicated 
by coincident alignments of unrelated pairs of stars often mimicking true wide binaries. Ac¬ 
curate kinematic measurements and statistical models are required to separate true wide 
binaries from coincident pairings. 

Wide binaries are also useful astrophysical probes, ideal for calibrating and testing em¬ 
pirical stellar property estimators. Gyrochronology, the process of estimating a star’s age 
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from its rotation period, is one of the most valuable techniques in stellar 


2007 h This can be used to derive ages for exoplanet host sta rs (e.g. 


2OI3I) or ultracool companions to intermediate mass stars (e.g. 


stellar astronomy 

(Barnes 

's (e.g. 

Walkowicz & Basri 

Duduv et al. 

2009) 

. These 


relations are often calibrated using open clusters of stars as anchor points. Each cluster 
provides a few hundred stars of different types, all of the same age. These discrete anchor 
points have the disadvantage that there are relatively few well st udied, older op e n clus ters 


with well determined rotation periods for their stellar members. 


Meibom et al. 


(1201511 re¬ 


cently provided rotation periods for the 2.5Gyr-old cluster NGC 6819. More held-age tests 
of gyro chronology relations will improve stellar age estimates even further, but each cluster 
still only represents a single age. Wide binary systems, however, cover a wide range of ages. 
While their age is not known independently as with clusters, they can still be used as a 
powerful test of gyrochonology at ages outside the range of known clusters. Wide binary 
systems are coeval systems and thus their components should have similar ages measured 
from gyrochronology. Kepler provides excellent photometry for rotation period calculation 
so any wide binaries of two Kepler target stars can be used to test gyrochronology relations. 
Gonversely, gyrochronology relations can be used to test if wide binary systems are true 
physical pairs. 

In this paper we present a proper motion survey of the Kepler held. We combine archival 
datasets from wide-held public surveys and target UKIRT observations with new Pan- 
STARRS 1 astrometry to produce accurate proper motions for Kepler target stars. This 
enables us to combine the excellent astrometric accuracy of a modern GGD-based sur¬ 
vey with the long time baseline of archival plate data. We then use our Pan-STARRS 1 
plus archival data proper motion catalogue, along with SED hts using Pan-STARRS 1 and 
2MASS photometry, to select a population of wide binary pairs where at least one compo¬ 
nent is a Kepler target star. We then use these binaries to test whether wide binarity has 
an ehect on exoplanet occurrence. 


2 DATASETS 

2.1 Pan-STARRS 1 data 


A 1.8 m high etendue survey telescope, Pan-STARRS 1 flKaiser et ahl 120021) recently com¬ 


pleted its full three and a half year survey operations on Haleakala on Maui in the Hawaiian 
Islands. This consisted of a suite of surveys of diherent cadences and depths. The Kepler 
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Wide binarity and planet occurrence in the Kepler field 5 


field was included in the Stt sur vey which covers the full sky north of 5 = —30° in five 
filters {gpi, rpi, ipi, zpi and ^Pi: lTonrv et ahll2012[) . Each filter has approximately six pairs 


of observations each separated by roughly half an hour over the course of the survey for 
any one point on the sky, therefore a typical object will have 30 pairs of observations. The 


data from Pan-STARRS 1 are astrometr i cally a nd photo r netric al 
using the processes o utlined in iMagnierl (120061) . iMagnierl (120071) . 


y reduced and calib rated 


(120121) and 


Schlaflv et ah 


Magnier et al.l ( 20131) . 


2.1.1 Astrometric accuracy 

There are three separate uncertainties in our Pan-STARRS 1 astrometric uncertainty model; 
centroiding errors caused by photon noise, the Pan-STARRS 1 internal systematic floor and 
the systematic floor caused by comparing Pan-STARRS 1 with other surveys. 

The astrometric centroiding error typically takes the form (Tpos = a? + Ku^ag where 
apos is our positional error in arcseconds, a^ag is the photometric error in magnitudes, a is 
our systematic floor and b is some conversion factor which depends on th e PSF shap e and 
whether the object lies in the source or background dominated noise regime (jKingill983l ). The 
Pan-STARRS 1 PSF takes a form which includes a free parameter (Magnier et ah, in prep.). 
Thus we cannot know b without knowing the exact PSF shape for every object. As this 
would be cumbersome and computationally impractical to implement, we take an empirical 
approach. We extracted all rpi detections in a 1x1 degree patch centred on R.A. = 18^42™' 
Dec. = 46°30™’ and measured the scatter in each object’s position about the mean position. 
We then plotted this as a function of the mean photometric error for the detections of that 
object (see Figure [I]). We found that a model of the form apos = y'^O.OlS^ -|- a^^g was a good 
£t to the data. Hence we adopt 15 milliarcseconds as our Pan-STARRS 1 internal systematic 
floor and set our photon-noise positional errors to be the same as the magnitude error (i.e. 
b = 1). 

To estimate the systematic error added by comparing Pan-STARRS 1 data with other 
surveys we again took another 1x1 degree test area (this time centred on R.A. = IQ^ld'" 
Dec. = 38°30™' as our first field has no SDSS data) and ex amined the positio nal scat- 


ter between Pan-STA RRS 1 and three external surveys (SDSS, 


Aihara et ah 


2011 


Skrntskie et ah 


2006 


; 2MASS, 


and P.I. UKIRT data of t he Kepler held processed using the same 


pipeline as UKIDSS data. 


Lawrence et ah 


20^ . Selecting only bright objects which lay 
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o 


mag 


a 


mag 


Figure 1. Astrometric scatter of rpi detections about mean object positions as a f unction of me an photometric error for the 
object’s rpi detections. The red line shows a least squares fit of the form dpos = \J^ the blue line shows our 

assumed astrometric model with a=15mas and 6=1. Green dots are the median scatter in each 0.01 mag bin. 


in the magnitude range of each of the surveys where the astrometric uncertainty is domi¬ 
nated by systematic errors (15 < r < 16 for SDSS and Pan-STARRSl, 11 < J < 15 for 
2MASS and UKIRT) we examined what the scatter about zero was for each pair of surveys. 
This provided us with six linear combinations of the systematic variances of the surveys. 
We solved these as a set of equations by singular value decomposition hnding that the sys¬ 
tematic variance of Pan-STARRS 1 data compared to external datasets is approximately 
38 milliarcseconds (including our internal systematic floor of 15 milliarcseconds). Once the 
internal systematic floor is subtracted in quadrature, we hnd a systematic error going from 
the Pan-STARRS 1 reference frame to some external reference of 35 milliarcseconds. Both 
external floor and internal floor are used later in our calculation of proper motions. 


2.2 UKIRT data 


Wemade use of observations which were made of the Kepler field using WFCAM flCasali et al 


200711 on the UK Infrared telescope which are now public. These data came from program 
U/09A/2 (PI Lucas) observed on July llth-13th 2009. Observations were done in the Jmko- 


band in a similar set-up the the UKIDSS Large Area Survey Jmkd observa 


crostepping, one telescope offset, 40 s total integration. 


Lawrence et al. 


ions (2x2 mi- 


2007) with reduced 


data catalogues available as FITS tables in the WFC AM Science Arch ive flHambly et al. 
20081) having been processed by the WFCAM pipeline fllrwin et al.l 120041) . 
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Wide binarity and planet occurrence in the Kepler field 7 


2.3 Public survey data 


We included astrometry from multiple pub 


the Two Micron A’ 
Survey (SPSS DR8 ; 


1 Sky Survey ('2M ASS; 


ic surveys in our wo rk. We drew data from 


Skrutskie et al. 


Aihara et al. 


2006 ). the Sloan Digital Sky 


flMouet et al 


201111 . the USNO -B digitisation of ph otographic plate data 


200311 and the WISE All-Sky Survey flWright et al 


201011 . The astrometry and 


photometry from these surveys were all extracted from the Vizier online archive service 
with the exception of WISE dataset which was downloaded from the NASA/IPAC Infrared 
Science Archive. To avoid spurious faint photographic plate detections affecting our proper 


motion hts we ignored any detections in USNO-B fa i nter t han Bj = 19, i?=19 anc 


This last cut is similar to that used in 


Deacon et al. 


( 20051) and IPeacon &: Hamblvl ( 2007 1 


In-7 S. 


For archival surveys we used the astrometric error estimates provided by the relevant 


survey archives with 


that found by 
uncertainty. 


he exc eption of SDSS. For this survey we used the error model used by 


Kraus fc Hillenbrandl (l2()07ll w hich features a 40 milliarcsecond systematic floor (similar to 


Pier et al. 


200211 and photon noise term that scales as the quoted magnitude 


3 CALCULATION OF PROPER MOTIONS 


We divided the Kepler held into 1x1 degrees chunks and in each of these areas extracted 
Pan-STARRS 1 data from the second large-scale reprocessing of the data (Py2l usi ng scripts 
written in the Desktop Virtual Observatory shell language (iMagnier et al.M2008h . We ex¬ 
tracted as our base catalogue the average catalog objects in the target area. This average 
object catalogue consists of photometric and astrometric properties derived from all obser¬ 
vations of each individual object and includes initial Pan-STARRS 1-only proper motions. 
As some objects will have proper motions that will take their archival detections outside our 
pairing radius of 1", we used the calculated Pan-STARRS 1-only proper motions to hnd the 
likely position of each object at the epochs of each archival non-Pan-STARRS observation. 
We then searched around these positions with a pairing radius of 1" to identify the objects’ 
detection in each survey. This procedure is vulnerable to spurious proper motions in the 
Pan-STARRS 1 database so we do not applying the proper motion “rewind” procedure for 
any object with a Pan-STARRS 1 proper motion > 10. 

To calculate our proper motions we constructed a covariance matrix taking into account 
that the Pan-STARRS 1 measurements will have three errors associated with them, photon 
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8 N.R. Deacon et al. 


noise, the systematic floor of relative Pan-STARRS 1 astrometry and the systematic error 
between Pan-STARRS 1 and other snrveys. The Pan-STARRS 1 measurements will have 
errors correlated with each other with a covariance of the square of the Pan-STARRS 1 
to other surveys systematic floor. No other survey data was given off-diagonal terms in 
the covariance matrix. We then calculated the proper motion for each object based on the 
measurements and our covariance matrix. 

We performed our fit in three stages. Having examined a number of preliminary fits 
we found that these were often offset by one or two Pan-STARRS data points which were 
significant outliers. Performing a 3a clip around our initial fit would be the most obvious 
solution. However this initial fit would be strongly affected by the outlier points and hence 
may exclude many valid datapoints which were outliers to this erroneous initial proper 
motion fit. Hence we performed an initial clipping removing any Pan-STARRS 1 data point 
which was more than a lOa outlier (i.e. a lOcr outlier in the quadrature sum of the number of 
standard deviations the object was an outlier from the fit in R.A. and Dec.). After this, we 
refitted our data and applied a more stringent 3a cut again only on the Pan-STARRS 1 data. 
This was followed by a final refitting of our proper motions. Figure [2] shows and example of 
the clipping and fitting process on one Aep/er target. 


3.1 Selecting the proper motions of Kepler targets 

Often there are multiple potential Pan-STARRS 1 matches for a particular Kepler target. We 
followed the following procedure to select the most appropriate matches. Firstly we selected 
only objects which had gpi, rpi or ipi magnitudes within 2 magnitudes of the g, r and i 
magnitudes of the Kepler target in the Kepler Input Catalog. This excludes matches with 
obviously spurious faint sources. We then matched each Kepler target with the predicted 
Epoch=1999.0 positions of all objects within 3" and implemented the following procedure 
for multiple matches. 

(i) We asked how many of those matches are within V and have r magnitudes in the KIC 
that agree with Pan-STARRS 1 rpi to within 0.2 mag. 

(a) If there is one such match we use this one 

(b) If there are multiple of these we use the one with the most measurements in its 
astrometric fit 
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Epoch 


Figure 2. An example of the fitting and clipping process for the star KIC 6878408. The initial fit is the red line, the second fit 
after the lOcr clip is the green line and the final fit after the 3(T clip is the blue line. The dashed lines represent the la confidence 
limits on the proper motion in each fit. Note points crossed out with red are excluded by the lOcr clip and those crossed out in 
green are excluded by the 3 a clip. 


(ii) If not we ask how many of the matches have r magnitndes in the KIC that agree with 
Pan-STARRS 1 rpi to within 0.2mags 

(a) If there is one snch match we nse this one 

(b) If there are mnltiple of these we nse the closest one 

(iii) Otherwise we nse the closest positional match 

3.2 Proper motion results 

Proper motions are presented in Table [U In onr remaining analysis we define significant 
proper motions as being more significant than 5a and good qnality proper motions as being 
calcnlated from 10 or more measnrements, having a USNO object within 6" and one PSl 
object within 3", being fainter than rpi = 14.5mag. (i.e. not satnrated), brighter than 
rpi = 19 mag and having a rednced l^ss than 4. This resnlts in proper motions which 
are not confnsed by crowding, satnration or ontliers and likely involving a time baseline of 
over 50 years. FignreOshows the distribntion of all objects in onr sample with nnsatnrated 
photometry (left) and photometry for nnsatnrated Kepler targets (right) in a rpi — ipi 
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Figure 3. The log^^Q density distribution of all unsaturated (rpi > 14.5 mag.) sources (left) and unsaturated Kepler targets 
(right) on a colour-colour plot. Note the concentration of Kepler targets around (0.2,0.5) in the Kepler targets plot. This is due 
to the large number of F and G stars in the Kepler target list. Note also the lack of an extension of the F, G, K locus after the 
stellar locus takes a turn at (0.5,1.2) with the onset of TiO absorption in M dwarfs. This suggests few Kepler targets in our 
sample are significantly reddened. 



Figure 4. The log^Q density distribution of reduced proper motion vs. colour for all sources with good quality 5cr proper 
motions (left) and targets meeting the same criteria. Note the clear white dwarf (marked WD), F, G, K dwarf (marked 

FGK), subdwarf (marked SD) and M dwarf (marked M) loci on the left-hand diagram. The right-hand plot is dominated by 
F, G and K stars, showing the known preference in Kepler target selection. 


Ys.gpi — rPI . This shows that there are relatively few Kepler targets which appear to be 
reddened beyond the stellar locns. Figure 0] shows reduced proper motion diagrams {Hr = 
Tpi + 5 + 51og^o/^) both all objects with good quality significant proper motions (left) 
and i^ep/er targets passing the same criteria (right). Clear F, G, K dwarf, M dwarf, white 
dwarf and subdwarf loci can be seen in the left-hand diagram, but the Kepler targets are 
disproportionally dominated by F, G, and K dwarfs. This is unsurprising as these were the 
primary target class in Kepler’’s search for Earth-like planets (iRatalha et al.ll2010l) . 
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Figure 5. The log^Q density distribution of proper motion error vs. magnitude (left) and reduced magnitude (right) 

for all sources. Note the discontinuity in proper motion error at rpi = 19. This is because we excluded matches with USNO 
detections at this magnitude to prevent matches with faint, spurious USNO detections. Note also both distributions become 
warped brighter than rpi = 14 due to saturation. Our reduced P^^t peaks at a value slightly below 1, indicating our 
astrometric errors are slightly over-estimated. 
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Table 1. Proper motions for unsaturated Kepler targets along with their Pan-STARRS 1 magnitudes. The penultimate column flags if this object passes the ’’good” proper motion bO 
test previously set out in the text and the final column shows Kepler Objects of Interest which have passed or failed our vial proper motion inspection. The number of objects within 3 
arcseconds and 6 arcseconds in both Pan-STARRS 1 and USNO are also quoted. Objects marked with f have multiple PSl matches within 3'' but the astrometric solution selected is . 
for an object within one arcsecond with the same magnitude as the Kepler target and is thus believed to be reliable. The full table is available electronically. ^ 




KID 

R.A. 

Dec. 

COS 5 


^meas 

9P1 

rpi 

ipi 

zpi 

ypi 

xl 

nSpsi 

n6psi 

'Ti^USNO 

§^6[/SiVO 

flag 


(Eq=J2000 Ep=2013.0) 

(mas/yr) 

(mas/yr) 

(mas/yr) 

(mas/yr) 


(mag.) 

(mag.) 

(mag.) 

(mag.) 

(mag.) 



Ct) 

o 



757450 

19 24 33.03 

+36 34 38.6 

12.08+1.65 

2.85+1.54 

59 

15.773+0.001 

15.087+0.001 

14.844+0.001 

14.737+0.001 

14.679+0.001 

0.485 

1 

1 

1 

s 

1 

1 

891916 

19 23 49.98 

+36 41 11.8 

8.79+1.61 

2.26+1.5 

57 

15.179+0.001 

14.691+0.001 

14.512+0.001 

14.439+0.001 

14.392+0.001 

0.626 

1 

2 

1 

a 

1 

1 

892718 

19 24 34.01 

+36 38 53.9 

-4.83+1.39 

-4.04+1.27 

65 

16.515+0.002 

15.842+0.001 

15.586+0.001 

15.451+0.001 

15.373+0.002 

0.553 

1 

1 

1 

2 

0 

892772 

19 24 36.80 

+36 40 43.8 

-8.16+1.3 

0.95+1.19 

74 

15.77+0.001 

15.05+0.001 

14.772+0.001 

14.627+0.001 

14.525+0.001 

0.72 

1 

2 

1 

1 

1 

1 

892832 

19 24 39.15 

+36 40 27.5 

-13.11+1.32 

-12.47+1.2 

81 

16.304+0.002 

15.624+0.001 

15.367+0.001 

15.252+0.001 

15.177+0.001 

0.326 

1 

1 

1 


2 

0 

892834 

19 24 39.20 

+36 37 39.3 

2.38+1.39 

2.28+1.28 

69 

15.919+0.001 

15.134+0.001 

14.863+0.001 

14.731+0.001 

14.657+0.001 

0.617 

1 

2 

1 


1 

1 

892882 

19 24 41.58 

+36 41 54.1 

1.89+2.47 

8.43+2.23 

23 

15.53+0.001 

14.862+0.001 

14.638+0.001 

14.542+0.001 

14.483+0.001 

1.363 

1 

3 

1 


1 

1 

892911 

19 24 43.13 

+36 40 14.1 

10.57+1.93 

0.18+1.86 

18 

16.081+0.001 

15.601+0.001 

15.433+0.001 

15.366+0.001 

15.314+0.002 

1.046 

2 

2 

1 


1 

0 

892946 

19 24 45.03 

+36 40 21.0 

8.44+1.85 

3.73+1.76 

25 

16.237+0.001 

15.804+0.001 

15.663+0.001 

15.622+0.002 

15.586+0.002 

1.285 

1 

2 

2 


2 

0 

893033 

19 24 49.82 

+36 40 03.8 

19.18+1.32 

-8.03+1.21 

77 

15.976+0.001 

15.197+0.001 

14.899+0.002 

14.776+0.001 

14.699+0.001 

0.442 

1 

1 

1 


2 

0 
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Figure 6. Comparison of our Pan-STARRS 1 proper motions with those from the UCAC4 catalogue ijZacharias et al]l2Q13h . 
Only points which pass our proper motion quality cuts are included. Red points show 5cr disagreements between our proper 
motions and UCAC4 in R.A (left panel) or Dec. (right panel). The majority of data points fall along a one-to-one relation 
between the two proper motions. A number of objects have low proper motions in our calculations but high UCAC proper 
motions 


3.3 Verification of proper motions 


To verify our proper motions we compared our selectio n of reliable proper motions with the 
values from the UCAC4 survey fIZacharias et al.ll2013l: see Figure [H]). This dataset does not 
use Schmidt survey plates so it represents an independent dataset. We found that our values 
generally compared well with the UCAC values. However, there was a small number (3.3%) 
of objects for which we measure much lower proper motions compared to UCAC. A visual 
inspection of a subset of these objects suggests that these objects do not have signihcant 
proper motions and often have another star at a 6-10" distance suggesting that these high 
proper motions are an artefact of object confusion in UCAC. As an additional check, we 


me 

in 


ividually in s pecte d the astrometric solutions for 3800 Kepler Objects of Interest listed 


Rowe et ah 


fl2015h for which we had a proper motion solution. We then made a by-eye 
judgement on how reliable the proper motion solution was compared to the distribution 
of astrometric points. We found that 94% of our proper motions passed this test. Of the 
objects which failed visual inspection, approximately three quarters had proper motions 
which would not meet our reliability cuts (high reduced few measurements, etc.). This 
KOI visible checking is listed as a parameter in Table [H 
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3.3.1 White dwarf Kepler Targets 

Proper motion catalogues with high levels of contamination due to spurious proper motions 
often present reduced proper motion diagrams with white dwarf loci which are heavily 
contaminated by FGK stars with erroneously high proper motion measurements. Hence 
the population of the white dwarf locus is a useful test of the reliability of a proper motion 
catalogue. To identify a sample of Kepler target white dwarfs we selected unsaturated objects 
which had good, signihcant proper motions, were redder than gpi — ipi = —2 (to exclude 
objects with erroneously blue photometry), met the condition > 10 x {gpi — ipi) + 12 
and whose proper motion solutions passes a visual inspection. These objects are shown in 
Table [21 Most of the intrinsically brighter objects passing this cut (those with low reduced 
proper motion) are sdO or sdB stars based on literature classifications (see Table [2]). This cut 
also misses intrinsically bright white dwarfs which will be distant enough that they would 
not pass our proper motion significance cut. Twelve of our fifty candidate white dwarfs have 
been previously identified as white dwarfs, dwarf novae or cataclysmic variables. We do not 
see signihcant contamination from objects spectrally classihed as FGK stars. 


4 IDENTIFICATION OF LIKELY BINARY COMPANIONS 
4.1 Probabilistic binary selection 


Stars in wide binary systems are typically found due to their common proper motion. Iden- 
tihcation of such pairings is complicated by the presence of nearby stars, which by chance, 
have the same proper motion as each other despite having no physical connection. These 
coincident pairings are often removed with hard cuts or using by-eye j udgements. Two bet¬ 


ter de hned processes fo r bina ry star selection come from the work of iLepine fc Bongiorno 


O 

o 

) and 

Dhital et al. 

(2010) 


(120in! ) . The former empirically dehnes a region of proper motion and 


separation space where coincident pairings are unlikely. The latter method uses a model of 
the Galaxy to dehne the probability that an unrelated star would have a particular proper 
motion and position difference from the target primary star. Here we dehne a hybrid of 
these two systems which calculates the probability that an object is a true binary based 
on a estimate of the background star population, but unlike Dhital (2010), our coincident 
pairing estimates are purely empirical. However unlike the two previously mentioned meth¬ 
ods, we take the likely distribution of the wide binary population into account as well. 
Similar methods of separating objects from a background population have been used to 
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Table 2. ^Ostenseiiet ay_j [201Clh . ^Ostensen et alj (I 2 OIIIJ) . 
ll2012h. teoffmeistei (Il966l~l. lOstensen et al. 

1'- Scaringi et alj ll2013f) . ^ ^Feldmeier et alj (l201lh . 


Greis^t alj ll2014h. 


(l2011bt), : 

Reed et alJ ' 

201 

if). iGirven et alJ ll201lh. ^McNamara et alJ 

(l2011af. ® 

Howell et al, 


2013f. ^Bischoff-Kim & 0 stenserj 

2011 

), 

^^Baran et alj ll201lh. 

i; 

Gianninas et alj ll201lf).^^Hoa: et alj i 

200C 

), 


Kepler^Y) 

Other name 

Literature 

classification 

Reference 

3353239 

2MASS J19364633+3825268 

sdB 

2 

3527751 

2MASS J19033701+3836126 

sdB 

3 

3561157 




3629120 

SDSS J190410.70+384518.1 

WD cand 

4 

3663415 




3751235 




4357037 

SDSS J191719.16+392718.8 

WD cand 

4 

5077438 

TYC 3123-1545-1 


14 

5544772 



5 

5951261 



5 

6042560 




6212123 




6672883 




7346018 




7594781 




7659570 

V* V344 Lyr 

DN 

6 

7797992 




8075923 




8077281 

2MASS J18501686+4358284 

B 

1 

8210423 




8244398 




8395780 




8420780 

GALEX J192904.6+444708 

DB 

7 

8490027 


CV cand 

8 

8626021 

WD J1929+4447 

DB 

9 

9071514 


CV 

10 

9139775 


DA 1 


9228724 




9391127 




9535405 


DA 

11 

9569458 


sdB 

1 

9818160 




10066680 




10081214 




10149875 




10290697 




10462707 


sdB 

2 

10590313 


LSPM J1910+4750 


11176123 




11191648 




11351218 



5 

11357853 


sdOB 

1 

11402999 




11509531 




11558725 

2MASS J19263411+4930296 

sdB 

12 

11604781 

2MASS J19140898+4936410 

DA 

2 

11822535 

WD 1942+499 

DA1.4 

13 

11911480 


DA/ZZ Ceti 

14 

12021724 

2MASS J19441275+5029393 

sdB 

2 

12353867 





i denti fy members of open c 


( 120041 ) ■ 


Kraus fc Hillenbrand 


usters usin gs their p r oper motion such as 
(120071) and 


Deacon &: Hamblv 


tion functions to the background population, we use the 


Sanders (1971). While these methoc 


Lepine fc Bongiornol (l2007l ) method 


s £t distribu- 


of offsetting the position of the target primary and then calculating the density of objects 
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16 N.R. Deacon et al. 


with similar proper motions, positions, distances and masses to the candidate secondary. 
Appendix ontlines the mathematics nsed in onr probability calcnlations. 


4.2 Testing the binary probability algorithm 

To test that onr method for binary association probability is accnrate we used a sam¬ 
ple of bright stars from our Kepler proper moti on catalogue. Thes e had proper motions 


from UCAC4 f Zacharias et al.l 12013 ) or PPMXL f Roeser et ah 2010) and photome 


tanc es and mass estima tes calculated using the methods of 
and 


Kraus et al. 


Kraus &: Hillenbrand 1 2(^7 ') 


fj2014bl ). We assumed a log-flat separation distribution flOpikI Il924j) Q 


flat mass ratio distribution flRaghavan et al 


ric dis- 


2010 l) and set the binary fraction wider than 


100 AU to be 25%. This latter number is similar to the number o 


limit found around nearby solar-type stars by iRaghavan et al.l (j2010|). We did not set a 


binaries wider than this 


lower proper motion limit and considered only pairs with observed separations below five 
arcminutes and projected separations below 10,000 AU. We used this technique to calculate 
companionship probabilities for these pairs and then flagged a subsample of these pairs for 
follow-up observations. 

A total of 26 wide binary ca ndidate systems (52 stars) were observed with the Tull Coude 


spectrograph flTull et al. 


19951) at the Harlan J. Smith 2.7m Telescope at McDonald Obser¬ 


vatory on the nights of August 8, 9, 10, and 11, 2014. Exposure times varied from 120s to 
1200s, which was sufficient to achieve a median SNR (per resolving element) of > 30. The 
resulting spectra cover 3800A to 10500A, with a gap from 4550A to 4750A due to internal 
reflection, and at a resolution of i? ~ 60, 000. Pair components were always observed within 
minutes of each other to reduce systematic sources of error (e.g., atmospheric changes, tele¬ 
scope flexure. Earths motion). Data were reduced using standard IRAF reduction tools, 
including bias and flat field correction, extraction of the 1 dimensional spectrum, and wave¬ 
length calibration. An additional wavelength correction was ap plied by cross cor r elatin g the 


observed telluric lines to a model atmosphere as described in 


Gnllikson et al. 


(I2ni4h . We 


cross-correlated the observations of our target pair components with each other to deter¬ 
mine the radial velocity difference. These values for our observed pairs (along with their 
binarity probability pun) are listed in Table |3l The errors in radial velocity difference come 


^ This is flatter than the typical log-normal separation distribution used for binaries closer than 100 AU but is appropriate as 
the binaries we identify in this work are typical a few thousand AU in separation 
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Figure 7. The observed radial velocity differences for a series of candidate binaries with a range of calculated binary proba- 
bilities. These objects had proper motions from UCAC4 and SED fits from the catalogue built for lKraus et alj ll2014bl') . Two 
of our pairs have a component which is a spectroscopic binary. This will induce offsets in their radial velocity. Note also the 
large fraction of objects which have high binarity probabilities and common radial velocities. 


from the cross correlation (with a Aa clipping) but do not include other potential sources of 
random and systematic error (e.g, additional undetected companions) and hence they are 
likely to be an underestimate of the error of the radial velocity difference. Figure [7| shows 
our results with two pairs flagged as containing a spectroscopic binary (i.e. are hierarchical 
triples). These were KIC 9085834 and KIC 11017620, which both produced a double-peaked 
correlation function or showed pairs of stellar lines. Nine of the pairs we classify as having 
Pbin > 0.6 have radial velocity differences less than 4km/s. One of these contains spectro¬ 
scopic binary while another three objects in the same probability range have radial velocities 
which are discrepant (excluding spectroscopic binaries). Of the hve low probability objects 
we observed all but one have discrepant radial velocities. This suggests that our method can 
reliably select true co-moving systems. 

4.3 Application to Pan-STARRS 1 Kepler targets 

The testing of our binary selection algorithm outlined in Section W?2\ only applied to objects 
with astrometry in UCAC4 and used only photographic plate and 2MASS photometry. We 
can widen our binary search by looking at fainter companions which have reliable Pan- 
STARRS 1 photometry. To do this we need a series of Spectral Energy Distribution (SED) 
templates using the process outlined in Appendix [Bl A by-product of this search allowed 
us to estimate accurate transformations from Pan-STARRS 1 colours to Kp magnitudes. 
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Table 3. Wide binary candidates with different binary probabilit ies Phin which had th eir RV differences me asured. Key: 
* Spectroscopic binary (SB2), fRapid rotator, Proper motions from ^Zacharias et all ll2013l') . feoeser et al. U 20 I 0 I) . 


KIC 

Ho, cos 5 

M5 

SpT 

Distance modulus 

r 

ARV 

Vbin 


(mas/yr) 

(mas/yr) 

(subtypes) 

(mag) 

(arcsec) 

(km/s) 


2696938 

5.2±2.0^ 

-18.3±2.0 

F6.7±1.5 

6.31±0.22 

5.4 

-0.55±0.85 

1.00 

2696944 

1.Oil.s'- 

-20.6±2.4 

F4.9±1.4 

6.53±0.22 




2992956 

16.5il.4'' 

50.2±1.4 

G6.5ii=1.45 

5.13±0.2 

16.9 

0.47±0.07 

1.00 

2992960 

14.5il.4^ 

51.8±1.4 

G6.9±1.45 

5.3±0.19 




4243796 

-1.2i0.7'- 

-18.7ihl.l 

K4.1ili0.6 

3.91±0.07 

154.3 

0.52±0.05 

0.00 

4346953 

-9.5i2.0'- 

-17.7ih2.5 

K7.2±0.15 

4.19±0.07 




5790787 

4.5i0.9'- 

10.6±0.8 

F6.9±1.45 

7.45±0.2 

27.6 

40.07±1.21 

0.70 

5790807 

s.iio.e' 

9.9±1.0 

F3.6±1.35 

6.52±0.23 




6934317 

-e.sio.s' 

-15.5±0.9 

G9.9±1.2 

6.46±0.12 

44.7 

38.29±0.11 

0.00 

7019341 

-il.Tio.e' 

-13.1±0.6 

G2.5±0.75 

6.66±0.05 




7090649 

16.Til.8^ 

13.6ii=1.8 

G7.0±2.6 

5.97±0.28 

9.3 

-2.11±1.91 

1.00 

7090654 

17.5il.2' 

13.3±1.2 

F7.2±1.85 

5.80±0.23 




7748234 

-6.2il.6'' 

3.9±1.6 

F8.2±0.2 

7.14±0.03 

38.0 

2.6±12.32 

0.74 

7748238 

-7.6i0.9' 

0.7±0.7 

F0.0±1.65 

6.8±0.17 




8123664 

-5.Oil.s'' 

-6.4±1.5 

F1.7±5.15 

5.96±2.5 

5.6 

8.63±0.62 

0.00 

8123668t 

-e.eii.s" 

-16.9ii=1.5 

F7.5ii=4.75 

7.70±0.49 




8619322 

-4.7i0.9'- 

-4.8±1.0 

K1.8±0.85 

3.62±0.11 

110.8 

-40.63±0.07 

0.51 

8683779 

-2.5il.3'- 

-6.3±0.9 

K4.4±0.6 

3.87±0.06 




9085833t 

-3.8i0.6' 

-15.9±0.8 

F4.8±1.35 

7.53±0.21 

20.1 

-0.22±5.35 

0.68 

9085834* 

-4.3il.6' 

-13.3ii=3.4 

G2.3±3.6 

8.36±0.32 




9579191 

-13.1i2.l' 

10.2±2.0 

G5.0±1.55 

7.05±0.15 

20.3 

-1.27±0.32 

0.0.86 

9579208 

-11.2i0.8' 

8.8±1.0 

F6.4±1.45 

6.39±0.21 




9595822 

-2.3il.O' 

-5.3±0.9 

K3.5±0.25 

4.67±0.03 

69.0 

29.64±0.11 

0.59 

9534041 

-S.Sil.l' 

-6.9±1.0 

K1.3±0.9 

4.92±0.12 




9655101 

-4.4i0.6'- 

-9.3±0.9 

K2.2±0.65 

4.73±0.08 

65.0 

2.7±0.1 

0.98 

9655167 

-4.7i0.7'- 

-8.7±1.1 

K3.0±0.45 

4.68±0.05 




9777293 

-o.sio.s' 

11.5ih0.7 

K1.6±0.85 

4.77±0.11 

48.9 

2.13±0.08 

1.00 

9777355 

-9.4il.3'- 

11.6ih0.8 

K2.0ili0.7 

4.79±0.09 




9912680 

11.6i2.0' 

1.0±0.6 

G0.4±2.7 

6.53±0.21 

7.1 

-0.5±0.27 

0.97 

9912690 

14.5il.3' 

2.5±1.2 

G0.9±4.6 

7.39±0.4 




10669568 

-0.1i2.0'- 

-10.9±3.4 

K4.2ih0.6 

5.68±0.06 

30.1 

60.48±0.17 

0.01 

10669590 

-10.6i3.5' 

-11.4±2.3 

K4.5±0.35 

5.91±0.04 




11017620* 

-9.4i2.l'- 

-19.0ii=0.5 

F7.1±1.95 

6.53±0.24 

8.9 

-18.93±0.6 

0.95 

11017626 

-4.8il.2'- 

-17.5±0.9 

F8.8±4.2 

7.21±0.38 




11551404 

iT.eii.s' 

12.6±1.1 

K2.6±0.6 

4.66±0.08 

55.7 

30.61±2.2 

0.00 

11551430 

11.2il.7' 

20.6ii=1.4 

K0.5±1.05 

4.88±0.13 




12253474 

-11.4i2.3' 

-8.1±0.5 

G0.7±4.65 

7.7±0.4 

18.0 

-1.02±0.17 

0.45 

12253481 

-5.4i2.2'- 

-7.0±1.0 

A6.6±5.05 

8.23±0.61 





This process is outlined in Appendix |P. We performed our SED fits using the same process 
as described by iKraus fc Hillenbrandl (120071 ) htting a reduced for all good data points, 


removing the most discrepant data point around the best fit (the £t with the lowest reduced 
X^) if it is more discrepant than 3a, rehtting, redoing our clipping of the most discrepant 
point with the same conditions as before and finally refitting. We do not clip if it would 
leave fewer than four good data points and we not do an initial £t if we have fewer than 
this number. Note that for the Pan-STARRS 1 photometry we exclude biters where the 
object is marked as extended and has fewer than four detections. We also exclude data from 
biters above the following saturation magnitude limits: gpi = 14.5, rpi = 14.5, ipi = 14.5, 
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Figure 8. A comparison between the measured spectral types and effective temperatures for our calibration stars and those 
derived from fitting. Note the objects shown here have at least four photometric measurement and exclude out M9.5 and LO 
calibration objects as we do not have an LO template to fit these against. In the upper plot a positive y-axis value indicates 
that our spectral type estimate from SED fitting is earlier than that derived from spectroscopy. Many of the M dwarfs have 
effective temperatures derived from their spectral types, hence the discrete temperatures in this region of the plot 


Zpi — 

used in 


4.5, yp^ = 14.5 an d Jmko = 10.5. This last saturation limit is the same as that 


Deacon et ah 


Kraus et al. 


(120091) . It should be stated that our SED templates (like those used in 
2014bl ) assume the target is a dwarf. Hence this method may miss or incorrectly 


match binaries containing one or more evolved stars. To test the accuracy of our SED 
htting method, we ran our calibration sample SEDs (as dehned in Appendix [B]) through our 
htting procedure. The results of these hts are shown in Figure HI This plot only includes 
astronomical objects with actual Pan-STARRS 1 measurement used for calibration not the 
synthetic Pan-STARRS 1 magnitudes used to model hot stars. 

The method produces reliable results for a range of spectral types, but with larger scatter 
for F and G type stars than K and M. Using a robust estimate of the scatter based on median 
absolute deviation we estimate that the standard deviation for our method is 2.5 subclasses 
for F and G stars and 0.6 subcl a sses for K and M stars. This equates to a scatter of 190 K and 

fj2ni4bl ) note that their SED htting is accurate to less than 


75 K respectively. 


Kraus et al 


one subclass. There is also the fact that the underlying calibration objects have measurement 
errors on their effective temperatures and/or those temperatures are derived from spectral 
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types which are typically quoted in discrete steps of half a subclass. For example our F and 
G calibration stars which are mostly drawn from iHuber et ahl (120141) have typical errors on 


their measured effective temperatures of about 110 K equating to 1.5 subclasses. This means 
that our quoted scatter in our measurements is likely to be around 180 K or 2.0 subclasses 
for F and G stars. We choose not to alter our scatter on K and M stars as we do not have a 
formal error on spectral type (and hence effective temperature) for many of these calibration 
objects. Our SED fitting procedure produces error estimates based on the calculated reduced 
distribution. However these are often extremely small confidence ranges compared to the 
scatter we have derived here. Hence when an object has a confidence range in spectral 
type that falls below our calculated empirical scatter, we substitute our empirical scatter 
measurements to estimate one sigma conhdence regions for our calculated distance moduli 
and bolometric magnitude. This results in a typical uncertainty in bolometric magnitude of 
0.4 mag. 


We included astrometry and SED hts for stars brighter than our Pan-STARRS 1 sat- 


uration limit of rpi < 14.5 by including astro metry from either 


Zacharias et al. 

2013 

) or the PPMXL catalogue ( 

Roeser et al. 

2010 ) 


he UGAG4 catalogue 


entry in both catalogues we selected the data source with the lowest quoted error on proper 
motion. To avoid losing objects around our saturation boundary we did not set any limits 
on magnitude, we simply excluded UGAG4 or PPMXL measurements of any object with an 
unsaturated P an-STARRS proper motion measurement in our catalogue. We then added 

fl2014bl ) giving us proper motion, spectral type and distance 


SED hts from 


Kraus et al. 


information for most bright stars in the Kepler held. 


We used our combined catalogue as a basis for our wide binary selection. From it we 
selected possible companions to Kepler targets separated by less than hve arcminutes from 
their primary and with distance moduli which dihered by less than two magnitudes. We 
then ran our binary probability algorithm using our quoted errors on distance modulus and 
proper motion for all our possible pairings. A histogram of our membership probabilities is 
shown in Figure [9l We have a large number of low probability pairings and an excess of high 
probability pairings with phin > 0.8. 
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4-3.1 Exclusion of giants 


One underlying assumption in our SED fits is that all the stars for which we derive distance 


moduli are c 


Mann et ah 


warfs . This is clearly untrue as many stars in the Kepler field will be giants. 


( 20121) showed that bright, red targets in the Kepler field are typically giants. To 


remove giants from our binary sam 

ole we set two cuts, one using o, r and Dr\^ photometrv 

from the Kepler Input Cataloj 

Both these use the separation 

me ( 

Brown e 

ah 

2011 

) and one using 2MASS photometry. 

Mann et ah 

(201^ 

!) showed between the late-type giant and 


dwarf populations. We define giants as either having J — Kg > 1 or having g — D^i > 0.25 
and g — r < 1.3. The hrst cut used here would exclude very late-type companion to a 
Kepler target so we do not impose this 2MASS related cut on our seconda ry stars. As an 
additional check we looked at the Hipparcos parallaxes (Ivan Leeuweru 120071) of stars in our 
binary sample. Of the 9 stars with measured Hipparcos parallaxes, none had distances which 
were significantly discrepant from our photometric distance estimates. 

We have endeavoured to minimise giant contamination but as the distance moduli for 
these binaries use the assumption that they are dwarfs, anyone using individual binaries 
from this sample should be aware that giant contamination may not have been completely 
eradicated and some giants may remain in the sample. We urge particular caution for the 
systems with primaries KIC 6367993 and KIC 10592818 as these have distance moduli less 
than one. 


5 THE WIDE BINARY POPULATION IN THE KEPLER FIELD 

After excluding giant stars and removing binaries closer than 6'!^ we were left with Kepler 
targets in 401 binary systems which have binary probabilities above 80%. We list these is 
Table m and show their projected physical and angular separations in Figure fTOl 


5.1 The ages of wide binary components 


Table m also contains rotation period estimates from lMcQuillan et al.l (120141) for stars where 
this was measured. For pairs with two Kepler Targets with measured rotation periods 
and separations grea ter than &' we estimated the ag es of both components using the age- 


rotation relations of 


Mamajek fc Hillenbrand 


(120081) . Rather than using measured B — V 


^ This cut was to remove binaries where the presence of a close companion star would cause errors in the determination of 
exoplanet parameters and rotation periods. 
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values (which do not exist for most Kepler targets) for each target we converted our mea- 
sured spectral types from SED fitting to colours using the field dwarf relation derived by 


Pecaut fc Mamaield fj2013l) . Figure ITT] shows the ages derived for our pairs. Most have rea¬ 


sonable agreement between their ages. Four of our 7 binaries agree within one a (roughly 
as expected) with one further binary having ages within two a. We would expect from our 
>80% probability binary selection criteria th at one or two of our seven pro posed binaries 
would be not true physical pairs. Additionally, iMamaiek fc Hillenbrandl (120081) find that four 
of the solar-type binaries in their test sample of 17 systems have ages which are discrepant by 
more than 0.3 dex. Hence our gyrochronology results do not indicate a larger than expected 
contamination in our binary pairs sample. 
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Table 4. Our candidate binary systems, the full table will be available electronically. Rotations periods from lMcQuillan et alJ ll201^) Key: ^ prope r 
^SED fit calculated using the methods of iKraus et alj proper motion/SED fit from this work, ^ proper motion from iRoeser et ^ 1120101) . 


motion from IZacharias et al. 


KIC 

KOI 

R.A. 

Dec. 

cos 5 


Distance modulus 

SpT 

Prot 

r 


Pbin 

WDS 



(Ep=2013.0, Eq=J2000) 

(mas/yr) 
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Pbin 

Figure 9. The distribution of binary probabilities for wide binaries with at least one Kepler target star. The arrow indicates 
the very large number of low probability pairings. The dashed line shows our 80% threshold for high probability binaries 


5.2 The planet-host fraction of wide binary stars 


To study how the presence of a wide binary companion affects the resulting planetary system 
we use our sample of wide binaries with companionship probabilities above 80% and which 
had passed our giant removal criteria. This sample of 401 multiples contained 529 Kepler 
targets. Next we crossmatched each Kepler targ ets in the sample w ith the most recent list 


of Kepler Objects of Interest (KOI) published by 


Rowe et al. 


(120151 ). These are divided into 


confirmed planets and planet candidates. If a binary system contained two Kepler target 
stars then we treated each Kepler target separately so a binary system where both Kepler 
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Figure 10. The properties of our binaries wider than 6" with binary probabilities greater than 0.8. The separations are shown 
in arcseconds (left) and AU at the calculated distance of the Kepler target (right). 
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Figure 11. Ages calculated from gyrochronology for binaries with two periodic Kepler targets as components. The lines join 
the two components of each pair with horizontal lines representi ng good agreement between the ages of both components. The 
error bars on the ages include the 0.05 dex errors suggested bv lMamaiek &: Hillenbrandl ll2008h . What appears to be a single 
data point is the binary pair KIC 10913758 and KIC 10913762 which our SED fitting classifies as having the same spectral 
type and which have near-identical rotation periods. 


targets were classified as a KOI we counted each of these as a separate planet-hosts. We 
restricted our sample to Kep/er targets with 10 < Kp < 16 and O.SMq < M < 1.5Mq as 
there are few Kepler target stars outside these ranges (where we assume a mass from our 
SED htting). Our separation greater than 6" cut meant we largely sample binaries with 
projected separations above 3,000 AU (see Figure ITOll . After our restrictions on mass and 
Kp magnitude we were left with 419 objects from which we identihed 10 KOI candidate 
planet hosts and 3 KOI conhrmed planet host (with two additional candidate planet hosts 
with masses outside our O.5M0 < M < I.SMq range). Hence 3.1±0.9% of our likely binary 
primaries are candidate or conhrmed planet hosts. No binary system where both components 
were KOIs was identihed. 

We also constructed a comparison sample of Kepler targets. These were 159152 objects 
which passed our giant exclusion criteria and had 10 < Kp < 16 and 0.5Mq < M < 
1.5Mq. We found that our binary primaries were generally brighter and lower mass than the 
comparison sample of Kepler targets in general. This is likely due to our method selecting 


^ We use the NExtSci Disposition for the candidate and confirmed planet host classifications 
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Figure 12. The properties of confirmed/candidate exoplanets around our binary components compared to a similar sample of 
Kepler targets. Our binaries exclude systems with separations below 


low-mass, nearby, high proper motion binaries preferentially as the number of coincident 
pairings with low-mass, nearby, high proper motion stars will be very low. To remedy this, 
we constructed 2D histograms for both our binary primaries and our comparison stars with 
bin sizes of Imag in Kp and O. 2 M 0 in mass. We then used these to construct a weighting 
factor for each comparison star based on its Kp magnitude and mass, 

Kbini^Kpy 
Kcompi^Kp^ Tn) 

where W is our weighting factor, Ncomp is the number of comparison stars in that bin and 
Nbin is the number of binary primaries in that bin. The median bin in Wm contained 12 
objects while the median bin in Ncomp contained 1804 objects. To estimate the number of 
planet hosts in our sample we crossmatched our comparison sample with the list of KOIs and 
then summed the weighting factors to calculate how many planet hosts would be expected in 
a sample of comparison stars of identical size, magnitude distribution and mass distribution. 
This led us to a comparison KOI rate of 2.4±0.5%. Hence we hnd no evidence for suppression 
or enhancement of the number of wide binary components which host exoplanets detectable 
by Kepler. We found that 1.0±0.5% of our binaries had more than one Kepler conhrmed or 
candidate planet while 0.5±0.1% of our weighted sample did. This is again a statistically 
insignihcant difference. We also examined the distributions of orbital period and planetary 
radius (see Figure [T^]) hnding no statistically signihcant differences between the candidate 
and conhrmed exoplanets around our binary components and around a similar sample of 
Kepler targets. 
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Our work finds no statistically significant difference between the fraction of wide binary 
components which are identified as being Kepler Objects of Interest and a sample of Kepler 
targets of similar brightness and spectral type. While these KOIs are not all conhrmed planet 
hosts, the fact that there is no significant difference indicates that in the regime where both 
our binary and comparison sample have detect planets (P<300 days) there is no evidence 
that a wide stellar companion with a separation wider than 3,000 AU will have a large effect 
on planet occurrence. Note that we do not claim that our comparison sample only contains 
single stars as many will have unresolved companions or companions which are too faint. 
This sample relies on a small number of binaries containing 419 Kepler targets. The ongoing 
Kepler K2 mission will drastically increase the number of stars in wide binaries probed for 
transiting exoplanets. This will increase the accuracy of the measurement of any possible 
effect from wide binarity. Additionally K2 will focus more on M dwarfs, objects largely 
excluded from our analysis he re. Finally i t shou ld be noted that our study does not disprove 


the theoretical hypothesis of 


Kaibetah 


(120131 ) that wide binaries can affect the orbits of 


planets. This effect could still be extremely significant for planets in wider orbits similar to 
the gas giants in our own Solar System. 


6 CONCLUSIONS 

We have used Pan-STARRS 1 data and archival datasets to calculate accurate proper mo¬ 
tions for stars in the Kepler field. This catalogue will be available to the community online 
and will improve the characterisation of Kepler target stars (i.e. giant-dwarf separation, kine¬ 
matic population membership etc.). By combining the proper motions with PSl-based SEDs 
to estimate the distance, we selected a sample wide binaries. To formalise this, we developed 
a statistical method to calculate true association probabilities for wide binaries including 
empirical estimates of the frequency of pairings of unrelated held stars. Our method was 
tested for a sample of bright stars showing that our likely binaries had low differences in 
radial velocity while our unlikely pairings had discrepant radial velocities. After excluding 
giant stars, we hnd 401 multiples containing at least one Kepler target which have binary 
companion probabilities of 0.8 or more and separations between six arcseconds and hve 
arcminutes. We hnd no diherence in the rate at which our wide binaries (typically wider 
than 3,000 AU) are identihed as conhrmed or candidate exoplanet hosts than a comparison 
sample of Kepler targets of similar brightness and spectral type. An investigation of those 
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binaries where both components have archival rotation periods shows that our binaries are 
roughly coeval indicating that our method is selecting true binaries. 


REFERENCES 

Aihara H. et al., 2011, The Astrophysical Journal Supplement Series, 193, 29 
Aller K. M. et ah, 2013, The Astrophysical Journal, 773, 63 

Baran A. S. et ah, 2011, Monthly Notices of the Royal Astronomical Society, 414, 2871 
Barnes S. A., 2007, 58 

Bastien F. A., Stassun K. G., Pepper J., 2014, The Astrophysical Journal, 788, L9 
Batalha N. M. et ah, 2010, ApJ, 713, L109 

Bischoff-Kim A., 0 stensen R. H., 2011, Asteroseismology of the Kepler Field DBV White 
Dwarf. It is a hot one 

Borucki W. J. et ah, 2010, Science (New York, N.Y.), 327, 977 

Brown T. M., Latham D. W., Everett M. E., Esquerdo G. a., 2011, The Astronomical 
Journal, 142, 112 

Gasali M. et al., 2007, Astronomy and Astrophysics, 467, 777 

Gushing M. G., Rayner J. T., Vacca W. D., 2005, The Astrophysical Journal, 623, 1115 
Deacon N. R., Hambly N. G., 2004, Astronomy and Astrophysics, 416, 125 

Deacon N. R., Hambly N. G., 2007, Astronomy and Astrophysics, 468, 163 

Deacon N. R., Hambly N. G., Gooke J. A., 2005, Astronomy and Astrophysics, 435, 363 

Deacon N. R., Hambly N. G., King R. R., McGaughrean M. J., 2009, Monthly Notices of 

the Royal Astronomical Society, 394, 857 

Dhital S., West A. A., Stassun K. G., Bochanski J. J., 2010, The Astronomical Journal, 
139, 2566 

Dressing G. D., Gharbonneau D., 2015, 21 

Dupuy T. J., Liu M. G., 2009, The Astrophysical Journal, 704, 1519 

Dupuy T. J., Liu M. G., 2012, The Astrophysical Journal Supplement Series, 201, 19 

Dupuy T. J., Liu M. G., Ireland M. J., 2009, The Astrophysical Journal, 692, 729 

Faherty J. K., Burgasser A. J., Gruz K. L., Shara M. M., Walter F. M., Gelino G. R., 2009, 

The Astronomical Journal, 137, 1 
Feldmeier J. J. et ah, 2011, 51 

Fressin F. et al., 2013, The Astrophysical Journal, 766, 81 


© 2015 RAS, MNRAS OOO.ITH391 


Wide binarity and planet occurrence in the Kepler field 29 


Gianninas A., Bergeron P., Ruiz M. T., 2011, ApJ, 743, 138 

Girven J., Gansicke B. T., Steeghs D., Koester D., 2011, Monthly Notices of the Royal 
Astronomical Society, 417, 1210 

Greiss S., Gansicke B. T., Hermes J. J., Steeghs D., Koester D., Ramsay G., Barclay T., 
Townsley D. M., 2014, Monthly Notices of the Royal Astronomical Society, 438, 3086 
Greiss S. et ah, 2012, The Astronomical Journal, 144, 24 
Gullikson K., Dodson-Robinson S., Kraus A., 2014, 7 

Hambly N. et ah, 2001, Monthly Notices of the Royal Astronomical Society, 326, 1279 
Hambly N. G. et ah, 2008, Monthly Notices of the Royal Astronomical Society, 384, 637 
Harris R. J., Andrews S. M., Wilner D. J., Kraus A. L., 2012, The Astrophysical Journal, 
751, 115 

Hewett P. G., Warren S. J., Leggett S. K., Hodgkin S. T., 2006, Monthly Notices of the 
Royal Astronomical Society, 367, 454 
Hoffmeister G., 1966, Astronomische Nachrichten, 289, 139 

Hog E., Fabricius G., Makarov V. V., Urban S., Gorbin T., Wycoff G., Bastian U., Schwek- 
endiek P., 2000, Astronomy and Astrophysics, 30, 27 
Howard A. W. et ah, 2011, 22 

Howell S. B. et ah, 2013, The Astronomical Journal, 145, 109 

Huber D. et ah, 2014, The Astrophysical Journal Supplement Series, 211, 2 

Irwin M. J. et ah, 2004, SPIE, 5493, 411 

Kaib N. A., Raymond S. N., Duncan M., 2013, Nature, 493, 381 

Kaiser N. et ah, 2002, in Survey and Other Telescope Technologies and Discoveries. Pro¬ 
ceedings of the SPIE, Volume 4836, Tyson, J. Anthony; Wolff S., ed., pp. 154-164 
King 1. R., 1983, Publications of the Astronomical Society of the Pacific, 95, 163 
Kraus A. L., Hillenbrand L. A., 2007, Astronomical Journal, 136, 2340 
Kraus A. L., Ireland M. J., Gieza L. A., Hinkley S., Dupuy T. J., Bowler B. P., Liu M. G., 
2014a, The Astrophysical Journal, 781, 20 

Kraus A. L., Ireland M. J., Hillenbrand L. A., Martinache F., 2012, The Astrophysical 
Journal, 745, 19 

Kraus A. L., Shkolnik E. L., Allers K. N., Liu M. G., 2014b, The Astronomical Journal, 
147, 146 

Lawrence A. et ah, 2007, Monthly Notices of the Royal Astronomical Society, 379, 1599 
Leggett S. K., Saumon D., Burningham B., Gushing M. G., Marley M. S., Pinfield D. J., 


© 2015 RAS, MNRAS 000.[Tti39l 


30 N.R. Deacon et al. 


2010, The Astrophysical Journal, 720, 252 
Lepine S., Bongiorno B., 2007, The Astronomical Journal, 133, 889 
Lepine S., Gaidos E., 2011, The Astronomical Journal, 142, 138 
Lepine S., Shara M. M., 2005, The Astronomical Journal, 129, 1483 
Magnier E. A., 2006, The Advanced Maui Optical and Space Surveillance Technologies 
Conference, E50 

Magnier E. A., 2007, in The Future of Photometric, Spectrophotometric and Polarimetric 
Standardization, ASP Conference Series, C. Sterken. San Francisco, ed., Vol. 364, Astro¬ 
nomical Society of the Pacihc, Blankenberge, p. 153 
Magnier E. A., Liu M., Monet D. G., Chambers K. C., 2008, Proceedings of the Interna¬ 
tional Astronomical Union, 248, 553 

Magnier E. A. et ah, 2013, The Astrophysical Journal Supplement Series, 205, 20 
Mamajek E. E., Hillenbrand L. A., 2008, The Astrophysical Journal, 687, 1264 
Mann A. W., Gaidos E., Ansdell M., 2013a, The Astrophysical Journal, 779, 188 
Mann A. W., Gaidos E., Kraus A., Hilton E. J., 2013b, The Astrophysical Journal, 770, 43 
Mann A. W., Gaidos E., Lepine S., Hilton E. J., 2012, The Astrophysical Journal, 753, 90 
Martin E. L., Cabrera J., Martioli E., Solano E., Tata R., 2013 
McNamara B. J., Jackiewicz J., McKeever J., 2012, B-STAR VARIABLES 
McQuillan A., Mazeh T., Aigrain S., 2014, The Astrophysical Journal Supplement Series, 
211, 24 

Meibom S. r., Barnes S. A., Platais L, Gilliland R. L., Latham D. W., Mathieu R. D., 2015, 
Nature, 517, 589 

Monet D. G. et ah, 2003, The Astronomical Journal, 125, 984 
Morton T. D., Swift J., 2014, The Astrophysical Journal, 791, 10 
Muirhead P. S. et ah, 2014, The Astrophysical Journal Supplement Series, 213, 5 
Newton E. R., Charbonneau D., Irwin J., Mann A. W., 2015, Accepted to ApJ, astro 
Opik E., 1924, Publications de L’Observatoire Astronomique de I’Universite de Tartu, 25, 
1 

Ostensen R. H., Bloemen S., Vuckovic M., Aerts C., Oreiro R., Kinemuchi K., Still M., 
Koester D., 2011a, At Last-A V777 Her Pulsator in the Kepler Field 
Ostensen R. H. et ah, 2011b, Monthly Notices of the Royal Astronomical Society, 414, 2860 
Ostensen R. H. et ah, 2010, Monthly Notices of the Royal Astronomical Society, 409, 1470 
Pascucci L, Apai D., Hardegree-Ullman E. E., Kim J. S., Meyer M. R., Bouwman J., 2007, 


© 2015 RAS, MNRAS OOO.ITH391 


Wide binarity and planet occurrence in the Kepler field 31 


ApJ, 673, 477 

Pecaut M. J., Mamajek E. E., 2013, The Astrophysical Journal Supplement Series, 208, 9 
Pickles A. J., 1998, Publications of the Astronomical Society of the Pacihc, 110, 863 
Pier J. R., Munn J. A., Hindsley R. B., Hennessy G. S., Kent S. M., Lupton R. H., Ivezic 
Z., 2002, 47 

Raghavan D., Henry T. J., Mason B. D., Subasavage J. P., Jao W., Beaulieu T. D., Hambly 
N. C., 2006, The Astrophysical Journal, 646, 523 
Raghavan D. et ah, 2010, The Astrophysical Journal Supplement Series, 190, 1 
Rayner J. T., Cushing M. C., Vacca W. D., 2009, The Astrophysical Journal Supplement 
Series, 185, 289 

Reed M. D. et ah, 2011, Monthly Notices of the Royal Astronomical Society, 414, 2885 
Roeser S., Demleitner M., Schilbach E., 2010, The Astronomical Journal, 139, 2440 
Rowe J. F. et ah, 2015, accepted to ApJS, astro 
Sanders W. L., 1971, Astronomy & Astrophysics, 14, 226 

Scaringi S., Groot P. J., Verbeek K., Greiss S., Knigge C., Kording E., 2013, Monthly 
Notices of the Royal Astronomical Society, 428, 2207 
Schlafly E. F. et ah, 2012, The Astrophysical Journal, 756, 158 
Skrutskie M. F. et ah, 2006, The Astronomical Journal, 131, 1163 

Tokunaga A. T., Simons D. A., Vacca W. D., 2002, Publications of the Astronomical Society 
of the Pacihc, 114, 180 

Tonry J. L. et ah, 2012, The Astrophysical Journal, 750, 99 

Tull R. G., MacQueen P. J., Sneden C., Lambert D. L., 1995, Publications of the Astro¬ 
nomical Society of the Pacihc, 107, 251 
van Leeuwen F., 2007, Astronomy and Astrophysics, 474, 653 

Walkowicz L. M., Basri G. S., 2013, Monthly Notices of the Royal Astronomical Society, 
436, 1884 

Wang J., Fischer D. A., Horch E. P., Huang X., 2015, The Astrophysical Journal, 799, 229 
Wang J., Fischer D. A., Xie J.-W., Giardi D. R., 2014, The Astrophysical Journal, 791, 111 
West A. A., Hawley S. L., Bochanski J. J., Govey K. R., Reid I. N., Dhital S., Hilton E. J., 
Masuda M., 2008, The Astronomical Journal, 135, 785 
Wright E. L. et ah, 2010, The Astronomical Journal, 140, 1868 

Zacharias N., Finch G. T., Girard T. M., Henden A., Bartlett J. L., Monet D. G., Zacharias 
M. I., 2013, The Astronomical Journal, 145, 44 


© 2015 RAS, MNRAS 000.[Tti39l 


32 N.R. Deacon et al. 


ACKNOWLEDGEMENTS 


The Pan-STARRSl Surveys (PSl) have been made possible through contributions of the 
Institute for Astronomy, the University of Hawaii, the Pan-STARRS Project Office, the 
Max-Planck Society and its participating institutes, the Max Planck Institute for Astron¬ 
omy, Heidelberg and the Max Planck Institute for Extraterrestrial Physics, Garching, The 
Johns Hopkins University, Durham University, the University of Edinburgh, Queen’s Uni¬ 
versity Belfast, the Harvard-Smithsonian Center for Astrophysics, the Las Cumbres Obser¬ 
vatory Global Telescope Network Incorporated, the National Central University of Taiwan, 
the Space Telescope Science Institute, the National Aeronautics and Space Administra¬ 
tion under Grant No. NNX08AR22G issued through the Planetary Science Division of the 
NASA Science Mission Directorate, the National Science Foundation under Grant No. AST- 
1238877, the University of Maryland, and Eotvos Lorand University (ELTE). We used public 
UK Infrared Telescope data tak en from program U /09A/2 (PI Lucas) this used the UKIRT 


Wide Field Camera (WFCAM; 


Casali et ah 


Hewe 


t et al 


(1200411 and 


200711 and a photometric system c 


(1200611. The pipelin e processing and science archive are described in 


escribed in 


Irwin et al 


Hamblv et al. 


(120081 1. The United Kingdom Infrared Telescope was operated by 


the Joint Astronomy Centre on behalf of the Science and Technology Facilities Council of the 
U.K. The authors wish to recognize and acknowledge the very significant cultural role and 
reverence that the summit of Manna Kea has always had within the indigenous Hawaiian 
community. We are most fortunate to have the opportunity to conduct observations from 
this mountain. The authors would like to thank Eddie Schlafly for discussions on correlated 
errors. 


APPENDIX A: DETERMINING COMPANIONSHIP PROBABILITIES 

In order to dehne a formal companionship probability, we dehne the following likelihood, 

0 = Ri ad, D, s, m) + 0/(/i, D, s, m) (Al) 

where (pc and (pj are the companion and coincident field distributions and are functions for 
the separation (s), proper motion difference (A/i), proper motion (/x) and distance modulus 
difference (AD) between any particular pair and the mass of the secondary component (m). 
The probability of a pair being real rather than a coincident pairing of unrelated held stars 

is 
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P = 


+ 


(A2) 


First let us examine the companion distribution, 


0c fbinfdet ^ 


-^-AD/2al- 


-g-ADV2(7|,- 
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Cs yoCL) 


(A3) 


Where A/i is the proper motion difference, is the quadrature sum of the proper motion 
errors, AD is the difference is distance moduli, ao is the quadrature sum of the distance 
modulus errors, s is the separation on the sky, m is the mass of the secondary and d is the 
calculated distance to the primary. The last two sets of brackets on the right each contain 
normalisation terms. We assume a power-law in physical separation but this method could 
be adapted to any underlying separation distribution or binary fraction. The separation 
distribution is normalised by Cg, 


Cs = 1/ f {sd) ^ ds (A4) 

J Si 

Where Si and S 2 are the inner and outer search radii of our survey. There are two free 
parameters to choose to complete this system, the hrst is the exponent x which dehnes the 
separation distribution. The other is the binary fraction f^in defined as the fraction of stars 
which have a companion in the region in which our separation distribution is defined. Next 
we must consider our observational selection effect caused by our choice of an inner and 
outer search separation. This is written in Equation I A3 1 as fdet and can be defined as, 

/*S2 C'f'max 

fdet = d {sdyds/ / (A5) 

d Si 

Where r is the physical separation in AU (dehned as sd) and our separation distribution 
is dehned in the region to Umax- Note the numerator of Equation IA5I is equal to the 
denominator of Equation IA41 Hence we can condense these two factors so that. 


/ T max 

r~^dr (A6) 

min 

The hnal term in Equation IA3I is the mass distribution of likely companions (/(m)) with 
the normalisation Cm- For this study we assume a hat mass ratio distribution. Hence as the 
mass of the secondary cannot be greater than the mass of the primary, the normalisation 
Cm becomes l/rup where nip is the mass of the primary (assumed to be the earlier-type star 
in the pair). 
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Next let us turn to the field distribution. The Galactic field star population is complex with 
the proper motion distribution particularly hard to p arameterise over a wide ran ge of proper 


Lepine &: Bongiornol fj2007() to measure 


motions. Hence we adapt the positional offset test of 
the likely contamination rate. Defining, 

(/)/ = Cf2TTs (A8) 

where cj is a the density of stars of similar absolute proper motion and has units of stars 
per (mas/yr.)^ per magnitude per sq. arcsecond per solar mass. The factor c/ is defined 
empirically by offsetting the pos itions of primary stars in o ur sample and repairing them to 
produce entirely confident pairs ( Lepine &: BongiornJ 2007 ). We do this and (for each of our 
actual pairs) measure the density of coincident pairs at the particular magnitude difference, 
proper motion difference and mass ratio of the actual pair. 


APPENDIX B: EMPIRICAL SEDS IN THE PAN-STARRS 1 
PHOTOMETRIC SYSTEM 

In order to calculate accurate photometric distances for stars in the Kepler field we require a 
series of empirical stellar Spectral Energy Distributions (SE Ds). There are multiple source s of 


empirical stellar SEDs with the most commonly used being iKrans fc Hillenbrandl (120071 ) for 


the SDSS/2MASS systems and IPecaut fc MamajekI ( 20131 ) for the Cousins/2MASS/WISE 


Tonrv et al. 

(2012 

) over the required range of spectral types (A-M). 

Duduv & Liu 

(2009) 


present synthetic colours for a range of spectral types b ased on init i al cha racterisation of the 


throughput of the Pan-STARRS 1 photometric system. 


Aller et al 


(120131) present a series of 


different SED templates for M dwarfs. These include templates for M and early L dwarfs in 
the star forming region Upper S co and early M dw arfs based on members of the Praesepe 


and Coma Her clusters. Finally 


Aller et al 


to be used templates for field stars. Note that 


(120131) presents indivi dual M9 and LO dwarfs 


Aller et al. 


(120131) used photometry from 


a previous processing version to our work so there may be subtle differences between the 
photometric systems. 

We began by taking the large collection of known dw a rfs in the Kepler field with spec¬ 
trally determined effective temperature from 


Huber et al. 


(I 2 OI 4 J ) and references therein. We 


converte d their effective temperatur es to spectral type using the dwarf effective temp e rature 


scale of 


Pecaut fc MamaiekI (120131) . For M dwarfs from the work of 


Martin et al. 


(1201311 
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we note that the effective temperatures in this work are derived from SED fitting so in¬ 


Pecant & Mamaiek 

(2013 

) scale. We supplement the obiects from the 

Huber et al. 

(2014) 


collection with late M dwarfs drawn from the compilation of 


Fahertv et al 


( 2009 ) as these 


a 2010.0 for 


West et al 


objects typically have high proper motions we corrected their positions to an epoc 
compa ring to Pan-STARRS 1 data. We also added M dwarfs from the sample of 
(|2008l) taking only bright < 16) and low extinction {Ag < 0.05) objects. We drew 

Pan-STARRS 1 photometry from the nightly science database and 2MASS data from the 
2MASS catalogue using pairing radii of 3 arcseconds. Objects were only incorporated into 
the fit in a particular PSl band if they had a magnitude error less than 0.05, more than 2 


We note that there are few objects in the 


(Schlafiv et al. 

201 

2 ). 


the 

Huber et a’ 

. (2 

014 

) catalogue hotter than 


6000 K which are not saturated in at least one Pan -STARRS 1 ba nd. Hence we supplemented 


our objects with synthetic colour s calculated fro m 


Pic 


Res 


(119981) and Dwarf Ar chives optical 
spectra and near-IR spectra from Cushing et al. ( 2005 ) and iRavner et ah ( 2009 ). The Pickles 


spectra overlapped with those from Rayner over the ypi band so we calculated mean optical 
colours relative to ypi in the AB system for each spectral type. We then calculated the 
colours relative to J 2 MASS by calculating ypi — J 2 MASS from the Rayner spectra, adding a 
term to take into account that ypi is in the AB system and J 2 MASS is in the Vega system. 

For the Pan-STARRS 1 bands we developed the following mechanism to calculate the 
magnitudes for each spectral type in each band. We began by fitting our gpi — ipi colours 
with a spline to produce a smooth distribution relative to spectral type. For spectral type we 
used ten subclasses per spectral classification.This is a departure from the traditional picture 
w here there are no K6 , K8 a nd K9 dwarfs. However we use the effective temperature scale 


of 


Pecant fc Mama jeld (120131) who include these missing spectral types in their temperature 


scale. Our gpi — ipi to spectral type fit is shown in Figure IBTI We then produced splines for 
relations between gpi — ipi and other colours; rpi —ipi, ipi — zpi, zpi — ypi and ypi — J2 MAss- 
These fits are shown in Figure IB^ We then used these colours at each gpi — ipi (and hence 
spectral type) to produce colours relative to the 2MASS J-band mag nitude. _ 


lese were then combined with the bolometric magnit udes from 


Kraus fc Hillenbrand 

(120071 ) and the bolometric correction and V — J colours of iPecaut fc MamajekI ( 2013 ). for 


objects later than M7 we used the spectral type - J-band absolute magnitude relations of 
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Figure Bl. The spectral type to gpi —ip\ transformation used to calibrate our SED templates. The black points are o bserved 
stars of kn own temperat u re/sp ectral type and the blue squares are synthetic photometry derived from the spectra of IPickl^ 
lll99d) and iRavner et ^ ll2009l') . The yellow line is a cubic spline fit to the data. 


Dupuv &: Lid (120121) . For near-infrared colonrs, we nsed — H-}.ma and Jqma ~ 


Ks 2 MAss colours derived from the bolometric corrections of 


Pecaut fc MamaiekI ( 2013 ). 


We also synthesised J2MASS - Jmko, H2MASS - Hmko and K.omasr - Kmko from the 


Rayner spectra to derive the MKO magnit udes flToknnaga 


et al. 


( 2012!) 


2002 ! ) of all our spectral 


spectral type - absolute 


types earlier than M7 and used the quoted iDupuv fc Liu 
magnitude relations for spectral types M7 and later. Our input sample included a handful 
of M9.5 and LO objects to anchor our hts at the cool end. We do not quote an LO template. 
Our hnal SED templates for each subtype between B9 and M9 are shown in Table | 
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9pi 'pi 9pi 'pi 


Figure B2. The transformation between Qpi — ipi and various other colours used to calibrate our SED templates. The black 
points are observed star s of k nown temperature/s pectral type and the blue squares are synthetic photometry derived from the 
spectra of |Pickl3 lll99d) and iRavner et ^ ll2009h . The yellow line is a cubic spline fit to the data. 
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Table Bl. The SED templates in the Pan-STARR S 1 and MKO photometric systems used in this work. The bolometric 
magnitudes and masses come from the compilation of iKraus Hillenbrandl J2007 ) and the effective temperatures are from the 

\eUt 


from 

Pecaut Sz Mamaiek 

1I2OI3I') 

|201J 

) andlOuDUv & Liu 

0 

s. 


dwar f effective temperatu r e scale is the dwarf scale from Pecau^^^Mamaj^ <l2013lh Our templates also use d the 2MASS colours 
Pecaut &6 Mamaiel^ J2^^) before M6 and from [Pupu^^^Liul 1120^1 . The reader should note that iPecaut &; MamaieS 
also include WISE colours. 


SpT 

9P1 

(mag.) 

’’PI 

(mag.) 

ipi 

(mag.) 

Zpi 

(mag.) 

ypi 

(mag.) 

Jmko 

(mag.) 

^MKO 

(mag.) 

^MKO 

(mag.) 

Mtol 

(mag.) 

K 

mass 

Mq 

B8.0V 

-0.31 

-0.03 

0.2 

0.34 

0.43 

0.02 

0.06 

0.16 

-1.0 

12500 

3.8 

B9.0V 

-0.05 

0.21 

0.42 

0.56 

0.64 

0.21 

0.22 

0.29 

-0.35 

10700 

3.35 

AO.OV 

0.34 

0.57 

0.77 

0.9 

0.97 

0.51 

0.51 

0.54 

0.3 

9700 

2.9 

Al.OV 

0.67 

0.87 

1.06 

1.18 

1.24 

0.78 

0.76 

0.79 

0.7 

9200 

2.65 

A2.0V 

1.02 

1.19 

1.36 

1.48 

1.54 

1.07 

1.03 

1.06 

1.1 

8840 

2.4 

A3.0V 

1.24 

1.38 

1.54 

1.65 

1.7 

1.21 

1.17 

1.2 

1.32 

8550 

2.267 

A4.0V 

1.41 

1.52 

1.67 

1.77 

1.81 

1.3 

1.24 

1.28 

1.53 

8270 

2.133 

A5.0V 

1.65 

1.73 

1.86 

1.96 

2.0 

1.46 

1.41 

1.44 

1.75 

8080 

2.0 

A6.0V 

1.87 

1.92 

2.03 

2.12 

2.15 

1.6 

1.54 

1.58 

1.92 

8000 

1.9 

A7.0V 

2.02 

2.03 

2.13 

2.21 

2.23 

1.67 

1.59 

1.63 

2.08 

7800 

1.8 

A8.0V 

2.2 

2.18 

2.26 

2.33 

2.36 

1.77 

1.67 

1.71 

2.26 

7500 

1.733 

A9.0V 

2.46 

2.41 

2.46 

2.53 

2.55 

1.93 

1.83 

1.86 

2.43 

7440 

1.667 

FO.OV 

2.67 

2.58 

2.62 

2.67 

2.69 

2.04 

1.92 

1.95 

2.61 

7200 

1.6 

Fl.OV 

2.84 

2.72 

2.73 

2.76 

2.78 

2.12 

1.97 

1.98 

2.75 

7030 

1.55 

F2.0V 

3.02 

2.86 

2.84 

2.86 

2.88 

2.2 

2.02 

2.02 

2.89 

6810 

1.5 

F3.0V 

3.35 

3.15 

3.11 

3.11 

3.14 

2.42 

2.23 

2.23 

3.13 

6720 

1.417 

F4.0V 

3.66 

3.42 

3.36 

3.35 

3.37 

2.63 

2.42 

2.43 

3.37 

6640 

1.333 

F5.0V 

3.96 

3.68 

3.59 

3.57 

3.59 

2.83 

2.6 

2.61 

3.61 

6510 

1.25 

F6.0V 

4.2 

3.89 

3.78 

3.74 

3.76 

2.96 

2.71 

2.72 

3.82 

6340 

1.223 

F7.0V 

4.48 

4.12 

3.99 

3.94 

3.96 

3.13 

2.88 

2.88 

4.03 

6240 

1.197 

F8.0V 

4.75 

4.36 

4.21 

4.15 

4.17 

3.31 

3.05 

3.06 

4.24 

6150 

1.17 

F9.0V 

4.91 

4.49 

4.33 

4.26 

4.27 

3.4 

3.13 

3.13 

4.35 

6040 

1.14 

GO.OV 

5.04 

4.59 

4.42 

4.34 

4.35 

3.47 

3.19 

3.18 

4.47 

5920 

1.11 

Gl.OV 

5.13 

4.66 

4.49 

4.4 

4.41 

3.52 

3.22 

3.21 

4.53 

5880 

1.085 

G2.0V 

5.16 

4.68 

4.5 

4.41 

4.42 

3.52 

3.19 

3.18 

4.6 

5770 

1.06 

G3.0V 

5.27 

4.78 

4.59 

4.5 

4.51 

3.61 

3.27 

3.26 

4.7 

5720 

1.053 

G4.0V 

5.38 

4.87 

4.69 

4.59 

4.6 

3.69 

3.35 

3.34 

4.79 

5680 

1.047 

G5.0V 

5.48 

4.97 

4.78 

4.68 

4.69 

3.78 

3.44 

3.42 

4.89 

5660 

1.04 

G6.0V 

5.62 

5.1 

4.91 

4.81 

4.81 

3.9 

3.54 

3.52 

5.03 

5590 

1.02 

G7.0V 

5.79 

5.25 

5.05 

4.95 

4.95 

4.03 

3.67 

3.65 

5.16 

5530 

1.0 

G8.0V 

5.94 

5.38 

5.17 

5.06 

5.06 

4.13 

3.76 

3.74 

5.3 

5490 

0.98 

G9.0V 

6.17 

5.58 

5.37 

5.25 

5.24 

4.3 

3.91 

3.89 

5.49 

5340 

0.94 

KO.OV 

6.39 

5.76 

5.54 

5.41 

5.4 

4.43 

4.03 

4.0 

5.69 

5280 

0.9 

Kl.OV 

6.67 

5.99 

5.75 

5.61 

5.59 

4.61 

4.18 

4.15 

5.89 

5170 

0.86 

K2.0V 

6.96 

6.21 

5.94 

5.78 

5.75 

4.75 

4.29 

4.26 

6.08 

5040 

0.82 

K3.0V 

7.32 

6.48 

6.17 

5.99 

5.94 

4.92 

4.41 

4.37 

6.32 

4840 

0.785 

K4.0V 

7.68 

6.75 

6.37 

6.16 

6.1 

5.04 

4.49 

4.43 

6.55 

4620 

0.75 

K5.0V 

8.09 

7.07 

6.62 

6.37 

6.29 

5.2 

4.63 

4.56 

6.68 

4450 

0.7 

K6.0V 

8.31 

7.22 

6.69 

6.42 

6.32 

5.22 

4.59 

4.49 

6.78 

4200 

0.665 

K7.0V 

8.51 

7.38 

6.8 

6.5 

6.4 

5.29 

4.64 

4.52 

6.89 

4050 

0.63 

K8.0V 

8.81 

7.66 

7.04 

6.74 

6.62 

5.5 

4.85 

4.72 

7.13 

3970 

0.617 

K9.0V 

9.11 

7.94 

7.29 

6.96 

6.84 

5.71 

5.07 

4.92 

7.36 

3880 

0.603 

MO.OV 

9.46 

8.27 

7.56 

7.21 

7.08 

5.93 

5.3 

5.15 

7.6 

3850 

0.59 

Ml.OV 

9.99 

8.76 

7.94 

7.55 

7.4 

6.24 

5.64 

5.45 

7.97 

3680 

0.54 

M2.0V 

10.8 

9.54 

8.54 

8.06 

7.88 

6.68 

6.07 

5.87 

8.44 

3550 

0.42 

M3.0V 

11.9 

10.61 

9.35 

8.76 

8.53 

7.28 

6.7 

6.48 

9.09 

3400 

0.29 

M4.0V 

13.27 

11.97 

10.41 

9.68 

9.39 

8.06 

7.53 

7.26 

9.92 

3200 

0.2 

M5.0V 

15.03 

13.69 

11.84 

10.97 

10.61 

9.19 

8.65 

8.35 

11.01 

3050 

0.15 

M6.0V 

16.58 

15.18 

13.03 

12.01 

11.56 

10.04 

9.46 

9.12 

12.06 

2800 

0.12 

M7.0V 

17.88 

16.38 

13.95 

12.73 

12.14 

10.56 

10.01 

9.57 

12.7 

2650 

0.11 

M8.0V 

19.05 

17.45 

14.73 

13.34 

12.62 

10.94 

10.36 

9.87 

13.13 

2570 

0.102 

M9.0V 

20.05 

18.33 

15.31 

13.83 

13.02 

11.21 

10.59 

10.09 

13.43 

2450 

0.088 


APPENDIX C: PAN-STARRS 1 TO KEPLER PHOTOMETRIC 
TRANSFORMATIONS 


Pan-STARRS 1 provides CCD photometry for both the original Kepler held and the ten 
K2 mission helds . Aep/er magnitudes are derived from a number of sources including Tycho 


Hog et al.l (120001), photographic plates and q . r, i photometry in the SDSS system from the 

fl2011ll. There are also additional 


Kepler Spectroscopic Classihcation Program 


Brown et ah 


datasets such as the Kepler INT Survey (KIS; 


Greiss et al 


20121) . In order to enable the 
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Table Cl. Color transformations from Pan-STARRS 1 filters to Kepler magnitudes. These take the form y = Cjo:* with 

the applicable colour range for each fit listed as Xrange 


ycolor 

^color 

Co 

Cl 

C2 

C3 

C4 

Grange 

Kp-gpl 

gpl-rpl 

0.0208986 

-0.871090 

-0.381325 

1.50797 

-1.12176 

-0.28-1.46 

Kp-rpl 

rp 1-ip 1 

0.0502748 

0.393252 

-0.650769 

0.126130 

-0.00725137 

-0.24-2.70 

Kp-ip 1 

ip l-2p 1 

0.131700 

2.48210 

-1.94587 

0.572845 

-0.0636664 

-0.15-3.13 

Kp - 2 p 1 

ip l-zp 1 

0.131700 

3.48210 

-1.94587 

0.572845 

-0.0636664 

-0.15-3.13 


calculation of Kepler magnitudes from Pan-STARRS 1 data we synthe s ised K epler magni¬ 
tudes for a series of stars using the dwarf spectral sequence of IPicklesI (119981 ) for types B9 
to M6. We supplemented these with seyen late M dwarfs from the Dwarf Archives spectral 
li brary □ with sufficien t spectral cov erage and thre e L an d T dwarfs from the compilation 


of 


Leggett et ah 


2010 1. Note for the 


Leggett et ah 


fj2010l) spectra we included objects with 


wavelength coverage until 4000A. This leaves the 3800-4000A region with no data in the 
Kepler filteT bandpass, however for objects this cool the emission in this region will be negli- 
gible. We synthesis ed colours using the Kepler bandpass @ and the Pan-STARRS 1 hlter set 


(ITonry et ah 


20121) . All colours were calculated in the AB system. Our colours are shown in 


Figure [CT] along with hfth degree polynomial hts. These coefficients of these hts are shown in 


Table ICll Figure ICll 


listed in 


Brown et ah 


so sh ows the synthesised Kepler magnitudes for the Kepler standards 


(120111 ) plotted against their observed Pan-STARRS 1 magnitudes. Ex¬ 


cept for two outliers which lie near the saturation point for Pan-STARRS 1, there is clearly 
fairly good agreement between our synthesised Aep/er-Pan-STARRS 1 and the Brown syn¬ 
thesised Kepler-real Pan-STARRS 1 colors. The small offsets are likely due to the differences 
in the ways the Kepler magnitudes were synthesised between the two sources. 


^ http://spider.ipac.caltech.edu/staff/davy/ARCHIVE/index.shtml 
^ http: //keplergo.arc.nasa.gov/keplerresponseyiresl.txt 
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Figure Cl. Synthesised colors fo r Pan-STARRS 1 to Kepler magnitudes for empirical spectra from l^ckle j lIlQQSl h Dwarf 
Archives and iLeggett et ^ ll2010l') (blue points). The red line shows th e fifth degree polyn omial fits listed in Table fUT] The 
aqua points are from the synthesised Kp magnitudes for the standards of iBrown et alj 1120111) and the observed Pan-STARRS 1 
magnitudes of these objects. 
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